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Abstract. Models with a scalar field coupled to the Gauss-Bonnet Lagrangian 
appear naturally from Kaluza-Klein compactifications of pure higher-dimensional 
gravity. We study linear, cosmological perturbations in the limits of weak coupling and 
slow-roll, and derive simple expressions for the main observable sub-horizon quantities: 
the anisotropic stress factor, the time-dependent gravitational constant, and the matter 
perturbation growth factor. Using present observational data, and assuming slow-roll 
for the dark energy field, we find that the fraction of energy density associated with the 
coupled Gauss-Bonnet term cannot exceed 15%. The bound should be treated with 
caution, as there are significant uncertainies in the data used to obtain it. Even so, 
it indicates that the future prospects for constraining the coupled Gauss-Bonnet term 
with cosmological observations are encouraging. 
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1. Introduction 

In recent years there has been a renewal of interest in scenarios that propose alternatives 
or corrections to Einstein's four dimensional gravity. These proposals are of differing 
origin as well as motivation: some are based on multidimensional theories in which 
gravity propagates in more than four dimensions (see for example pQ), others on scalar- 
tensor couplings, that could give cosmologically observable effects [2], or violate the 
equivalence principle [3]. Another class of theories modify gravity by adding terms in 
the Lagrangian that depend on quadratic combinations of the Riemann tensor (see for 
example jl]), such as R 2 and R^R^, yielding higher than second-order field equations. 
Among the applications of such a generalisation we notice that these terms have been 
used in the past to model inflation [5] and dark energy [6] . As it is well-known, there is 
a unique quadratic combination of the Riemann curvature tensor that, if added to the 
usual Einstein-Hilbert action, does not increase the differential order of the equations 
of motion [TJ E]- For historical reasons this is called the Gauss-Bonnet (GB) term as its 
geometric origin can be traced back to the Gauss-Bonnet theorem regarding the Euler 
characteristic of two dimensional surfaces [9]. 

In classical gravity theories involving dimensions higher than four there is no 
good reason to omit it (apart from complication). Studies of higher dimensional 
gravity, mostly in the context of braneworlds, have shown that such theories have 
many surprising properties. Without trying to be exhaustive, we can refer to black 
hole physics (e.g. [ID]), gravitational instabilities (e.g. [H]), and the dynamics of co- 
dimension one and higher braneworlds (e.g. [T2]). 

In four dimensions the Gauss-Bonnet combination reduces to a total divergence 
and, as such, is dynamically irrelevant. However, in the case of a scalar-tensor theory 
(see for example [13] for a detailed and general discussion), the Gauss-Bonnet term 
couples to the scalar sector and four dimensional gravity is modified (see also [H] for 
the case of Lorentz and Chern-Simons terms). A coupled Gauss-Bonnet term should 
therefore be included in the most general second-order scalar-tensor theory. Such a 
coupling is manifest if we consider, for example, the case of when the scalar field is a 
modulus field originating from a Kaluza-Klein (KK) toroidal (i.e. flat) compactification 
of a 4 + iV dimensional theory. This is explicitly shown in section [21 There it is also 
shown that in order to be consistent with the KK truncation, higher-order scalar field 
terms must also be included; a fact that is very often neglected in the literature. 

The above discussion about uniqueness of the Gauss-Bonnet term is, strictly 
speaking, valid only at the classical level. However, additional motivation stems from 
low-energy effective string actions, such as the one for heterotic string theory. At tree 
level (for string coupling) this includes the Gauss-Bonnet term as the leading order and 
unique (up to that order) ghost-free a' correction to the Einstein-Hilbert term [15] (see 
for example [16J for cosmology and braneworld models emanating from such actions). 

Since the coupled Gauss-Bonnet term is to be naturally included in a scalar-tensor 
theory, an important question arises: how can current observations constrain such 
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models? It has already been shown in the literature that at the background level a 
coupled Gauss-Bonnet term allows for a viable cosmology [17] . In this article we address 
the question of the evolution of linear cosmological perturbations, and see whether 
interesting constraints can be derived with present or future data. Another, and maybe 
more stringent way to set up constraints is to look for modifications of gravity on local 
scales; parameterized post-Newtonian (PPN) analysis [HI [19] requires a next to leading 
order study of slow-moving sources around the Earth's gravitational field, or of solar 
system bodies. Such a study would necessarily be dependent on the details of the Gauss- 
Bonnet coupling, since it is dimensionful, and as we will see would be independent of 
the cosmological expansion. In fact in this case the relevant background would be that 
of flat spacetime, or of a modified Schwarzschild metric. This study is beyond the scope 
of the present paper (see [20] for a first approach in this direction) and here we will stick 
to cosmological constraints which are as model free as possible. 

As a first and most conservative approach, and in order to disentangle the Gauss- 
Bonnet effects from the pure scalar sector of the theory, we place ourselves in the Einstein 
frama^]. Since we assume the scalar field is driving an accelerated expansion, we also 
endow the field with a sufficiently flat potential in a cosmological background. 

In section [2] we describe the class of theories we will be using, and give the 
cosmological field equations in suitable notation. We also derive the form of the theory 
when the scalar field is of Kaluza-Klein origin. In section [3j which is the backbone 
of this work, we write the cosmological perturbation equations. We then analyse the 
field equations for small Gauss-Bonnet coupling, and for a slowly rolling scalar field. 
Present-day constraints are discussed in section [U We conclude and discuss our results 
in section [5l 

2. Second-order scalar-tensor theory 

Let us consider the following Lagrangian 

lo7rG 

where 

£ (c) =ei(0)/: G i ? + 6(0)G^V^V> + e3(0)[V0] 2 V 2 + e4(0)[V0] 4 (2) 
£-gb = R 2 — ^Rfj,vR^ u + Rai3fj,vR al3 ^ v ■ (3) 

G^ v is the Einstein tensor and Cqb is the Gauss-Bonnet term. The coupling constant a 
has dimension mass -2 . Note that we place ourselves in the Einstein frame. Under this 
hypothesis the above expression (CQ) is the most general second-order ghost-free scalar- 
tensor Lagrangian (with uncoupled matter sector). In this article we will be interested 
in the observational effects coming purely from higher-order gravity, and so we will take 

| We will assume however that the energy- momentum tensor is conserved i.e., that the Einstein frame 
is the physical frame. 
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<5£ mat /<50 = 0. This hypothesis, along with the fact that we place ourselves in the 
Einstein frame, ensures that all deviations from standard gravity will originate from the 
£( c ) sector in 

The form of the functions 6 will be related to the origin of the scalar field. 
One possibility is the compactification of extra dimensions, as was mentioned in the 
introduction. Consider the general 4 + N dimensional Lagrangian of pure gravity, 

C 4+N oc (R - 2A + aC GB ) (4) 

which is of second order in the curvature operatojjjl and yields second-order field 
equations, which are divergence free and have well defined and stable perturbations 
around the vacuum. We use the metric anzatz, 



e^ ix) g^(x)dx»dx u + dX a dX a e 



-2[i4>(x)/N 



(5) 

where p 2 = 2N / (2 + N) and dX a dX a is the Euclidean metric, and also locally that of 
the iV-torus. Compactifying the N extra dimensions then gives the above theory ([I]), 
with £f 1 oc U oc e^, and 
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-X 



2 + N* 1 ' 1,0 v /iV(2 + iV)" 1 ' ^ (2 + N)N 
The scalar field plays the role of the overall size of the ./V-torus as we can see from (j5J) . 
We see that typically all of the four possible ghost-free quadratic-order gravity terms 
will be present, and not just the Gauss-Bonnet term as is often assumed (see also [2~T]). 

The general field equations have been given in [22J. Here we write down the 
equations for a Friedmann-Robertson- Walker (FRW) metric, including pressureless 
matter. We set the conformal Hubble function Ti = aH(t) and adopt the e-folding 
time r\ = log a as time variable: 



6 



2V2(2 - N) 



i • 



4(1 - N) 



6 



(6) 



3H 2 = 8nGpa 2 



Ua 2 + 



n + 1 



8ttG 



pa 



+ 2- 



2 

Ua 2 



n 



2i2 



aH 4 



m 2 m 2 

aH 212 



2a? 

~ 2a 2 



2^[-9^-6^ 2 -(3U-Q^ 3 



(7) 



+ 




12& - 36 

_ aH 2 '4> 
~ 2a 2 

12# - 960 - 9&« 

aH 2 



4ei - 26 



-6 



(8) 



36 + 660 + 2(36 



e 3 )0 2 (9) 



3^+ +(3e;-c 

(dots denote d/dr], primes denote d/d<f)). We now introduce the dimensionless variables 
%nGpa 2 2 Ua 2 „ 4H 2 

{>K — — , 




m 2 



2^ 6 



3ft 2 ' 

^ ISC 



fix 



ft. 



6a 2 



a0 4 (36-e 3 ), (10) 



§ For N > 2 to keep all generality [7] we would have to add higher-order curvature invariants. To 
understand their geometric origin see [5]. 
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which are subject to the Friedmann equation ([7j), so 

£a = i. (ii) 

Defining fi^ = 1 — f2 m , equation (jSj) gives the acceleration of the universe to be 
n/H = -[3Q^W0 + l]/2 where 

f^W<£ (l - + 2 ^ 2 + ^3 \ = _Q p _|_ ^ 






and gives the effective equation of state for all the ^-dependent terms, including 
the quadratic-order gravity ones. It is clear that the density fractions Qi quantify the 
deviation from standard gravity and standard cosmology. It is important to notice 
that since the & couple to TC, the densities fli will disappear from the field equations 
whenever 7i is zero; in other words for a(t) = constant. This signifies that we need 
a FRW background which is neither Minkowski nor de-Sitter in order for the higher- 
order terms to affect the cosmological evolution. This spells out why one needs to go 
beyond linear order in order to find local effects of the Gauss-Bonnet term, at least for 
a flat spacetime. With this in mind the goal of this article will be to express observable 
quantities in terms of the cosmological 



3. Cosmological perturbation equations and the linear post-Newtonian 
limit 

We derive in this section the cosmological perturbation equations; more precisely, we 
derive the limit for a weak gravitational field (i.e. linearised perturbation equations), 
with slow time- variation (we specify below how slow), and at small scales (with respect 
to the horizon's size). We refer to this as the linear post-Newtonian (LPN) limit. It is 
convenient to use the longitudinal gauge metric written directly in r\ = log a: 



ds 2 = e 2 " 



2 



;i + 2*A + (l-2$)dx f dz« 
ri 



(13) 



The field \1/ is the Newtonian potential and $ is the leading-order, spatial post- 
Newtonian correction which will permit the calculation of the stress-anisotropy. We 
always assume that the effective mass of the field is very small, which is what is 
expected for a dark energy field; the precise requirement is that the length scale 1/m^ 
be much larger than the typical scales of the experiment used to set the constraints. 

In practice, the LPN limit amounts to considering A$ 3> 7i 2 x ($, $, $) and 
similarly for the other gradient terms. For a plane wave perturbation of wavelength A 
we see that if 2 $ is much smaller than A$ when A 1/H. The requirement that $ is 
also negligible implies the condition 

d -^«(XHr\ (14) 
a log a 
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which clearly holds in perturbation theory for any reasonable perturbation growth as 
soon as A <C 1/H (notice for instance that the gravitational potential $ is in fact 
constant in a standard matter dominated universe, or in a spherical collapse, and slowly 
varying in most cosmological models, e.g. those with a cosmological constant). The 
same arguments apply for $ as well as the other metric potential and also for 5<f). 
The (7777) and (ij) components of the gravitational field equations, and the scalar field 
equation, can be written in the following compact fashion, 



-H 2 5tt m = AA$ + B 



Al 



= CA$ + D 



BA$ =DA$-E 



A50 



Ac 



(15) 
(16) 
(17) 



The background coefficients are given by 



.4 



C 



D 



E 














H 


. 2 + 2n 2 


1 - 


4 














1 - 



3 



fii + 2 



^20 

36 



1 + n + 



n 2 + s 



+ 



1 



3 6 

^0 
96 



+ AVL 4 



fi 3 + 4fi 4 .(18) 



For completeness we also include the (Oj) component of the gravitation field equations 
in the same limit 



-H 2 9Q m = (B - A) A^ + F 



Al 



AA&-B 



Al 



(19) 



2 

Finally, the energy-momentum conservation equations are as usual 5 = —9 and 
9 + 9 + (H/H.)9 = —A^/fTi. 2 , where 5 = Sp/p is the matter density contrast and 9 
is the divergence of the matter peculiar velocity field. 

In Einstein gravity, the perturbation equations in the LPN limit would give the 
usual expression of the Poisson equation (|T5l) and the condition \1/ = $, i.e. the vanishing 
of the stress anisotropy a = ^ — $ (1161) . To quantify the deviation from Einstein 
perturbations we will express our results by two observable quantities: the time variation 
of Newton's constant, G/G, and the fractional anisotropic stress 7 = Making use 

of the perturbation equations we find 
EA + BD 



7 



D 2 + CE 



(20) 
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Poisson's equation is 



Gr 2 



(21) 



where the variable Newton's "constant" is 



G*(r]) _ D 2 + CE 



(22) 



G ~ A 2 E + 2ABD — B 2 C ' 



From Poisson's equation, we see that the growth of the matter fluctuations is 
immediately given by the equation 



In the absence of any quadratic-order gravity terms (a = 0) one has 7 = 1 and 5<p = 0, 
returning a constant effective gravitational coupling (G* = G). Notice in passing that 
the corrections to gravity cancel in the limit 7i — > 0, which shows that the contribution 
from the coupled Gauss-Bonnet term vanishes at the linear level for this limit. 

Actual experiments set up to constrain the variation of G (and the value of 7) are 
generically local experiments set-up in the immediate vicinity of the local gravitational 
field of the Earth and the Sun [THJ HH]- Also cosmological constraints on G from 
primordial nucleosynthesis constrain the variation in the expansion rate induced by 
a variable G in an otherwise standard Friedman equation (see e.g. [23]). In the present 
case, however, the Friedman equation does not contain a variable G since we did not 
include the Brans-Dicke term. Had we done so it would be difficult to disentangle 
the effect of the scalar field from the Gauss-Bonnet correction. In the actual state of 
experiment it turns out that the only way to observe our modified G is to constrain the 
form of the Poisson equation from cosmological observations. 

The Poisson equation of the GB theory is similar to the one in Brans-Dicke gravity 
which (in the Einstein frame) arises from the coupling of the scalar field to matter or 
(in the Jordan frame) from scalar corrections to gravity as a generalised Lagrangian 
C(<p, R). It is clear that the specific form and time dependence of the correction to the 
Poisson equation G*(rj), depends on the form of the Lagrangian. We expect 7 — 1 and 
G* to be small when f2j are small, and so we will now examine several different limiting 
situations for which this is the case, either as a result of small a, or slow- roll (small 0). 

3.1. Small-a limit 

Let us begin with the small-a limit. If we include terms which are of linear order in a, 
we find that 8(f) and a are no longer zero. Solving the perturbation equations ffl6|) and 
( TT71) . we find that in this limit the stress anisotropy parameter is 




(23) 




(24) 
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In general, the constraint on the Gauss-Bonnet densities f2j will depend on the specific 
dynamics. Similarly, the gravitational constant is 

m m c_ l+ ^_t_^ +tk+ol ^ (25) 

There are also in principle O(aA<50) corrections to the expression for G*, although we 
see from the scalar field equation (fTTI) that A<50 = 0(a), and so its contribution to G* 
is of quadratic order in a. The higher-order gravity terms Q3 and Q4 only affect G* via 
50, and so similarly their contributions are quadratic order in a. This is despite the 
fact that their contributions to the Friedmann equation are of the same order as Q\. 

Differentiating the expression (|25j) . and using the field equations (in the small a 
limit) to simplify it, we find 

G* 




(26) 

Again, to derive a precise bound on f2 lj2 from this expression requires the dynamics 
of to be specified, i.e. a choice of potential and couplings. However, unless the couplings 
and the potential are unnaturally steep functions of 0, we can assume that all the ratios 
U"/U', and similar are 0(1); then the expression simplifies to 



G* 
G* 



h + b 2 (j) + 6 3 2 + 5^ \Q 1 + 64 + 6 5 + 2^ fi 2 + 0(a 2 ) , (27) 



where all the fej's are quantities of order unity. 

We see that the coefficients of Qi^ are smallest in absolute value in the slow-roll 
limit, when <f)/<fi, <fi 1 (assuming there is no "chance" cancellation of the coefficients). 
With this in mind, we will now turn to cosmologies where the only assumption is slow- 
roll for 6. 



3.2. Slow-roll limit 

An assumption which is justified in the context of dark energy models is to suppose 
that the scalar field is in a slow-rolling regime. This allows us to specify to some extent 
the background dynamics of the field. Slow-roll is realised for 0/0,0 1. Under this 
hypothesis we can solve the field equations keeping leading-order kinetic terms. Our 
aim is to generically obtain the role of the higher-order terms of the Lagrangian, in the 
slow-rolling regime, without any further assumptions. In particular we do not take a to 
be small, in contrast to the previous subsection. 

The Friedmann equation truncated to linear order gives simply 

1 = VL rn + VLp + ^ + O(0 2 ) . (28) 
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In the same way the acceleration field equation (TT2"]) reduces to 

^ = 1 - - + O(0 2 ) (29) 

where we have used (1281) to trade flp for Q m . Note that a negative Gauss-Bonnet density 
fix can assist late time acceleration, and a positive one can suppress it. Substituting 
and (I2"9"j) into the scalar field equation Q we obtain, 



2^(l-a n ) + (2-3ft m )ft 1 

(2 - ft m )0 2 + 2(2 - 30 m )0 2 - ^-(1 + a m )^! - 6fi m fi 2 



+ 0(0^,0/0) = 0. (30) 

The leading-order part of the scalar field equation (130]) constrains the validity of the 
slow-roll limit. It is self consistent if, and only if, 

U' _ (3n m - 2) gi 2<H 2 (30 n -2) 

[/ 2(1 -fi m ) a 2 (l-^ m ) ' 1 J 

This shows that for 0^0, U'/U remains finite. Hence, in contrast to conventional 
first-order gravity, it is possible to have slow- roll even when U is not particularly flat. 
We see that to next order, ( 130]) gives us Q 2 with respect to Q\ and the scalar field 
parameters, 



2 



4(1 - n m )n 2 = _ 2 ( 2 "m)(l fim) ^ + ^ 2 9Q m n 2 (33) 

Taking the LPN limit we obtain 

7 = i- , njp-auq-aru) + 0(#> M) (33) 

20 2 + 4fi 2 (l-fi m )+fi 2 (2-3fi m ) 2 



and 



G* . fi 2 (l-3fi 



2 

m y 



- . , + 0(0,0/0). (34) 

G 20 2 + 4fi 2 (l-fi m ) + 3fi 2 (l-2fi m ) 

One could further substitute Q 2 from (132]) but the resulting expressions are more 

cumbersome. The time variation of the gravitational coupling to the same order is 

then equal to 

k=- 6 w^wi 3a - +3 ^- 2 )i i -B- (35) 

In contrast to the small-a limit, this time the 50 perturbation does give significant 
contributions to 7 and Furthermore the corrections to 7 and are large even if 
we take — > 0. This means that taking a slow-roll solution is not sufficient to satisfy 
observational constraints. 

As with the small-a limit, the contributions from £ 3 and £4 do not feature in the 
expressions for 7 and G>, although this time it is due to the smallness of rather than 
50. Finally, we notice that a chance cancellation at first order of both non-Newtonian 
effects 7— 1 and G*/G (and G*/G—l) in (133]) and (|34j) occurs for Q m = 1/3, a remarkable 
value indeed. In this case it is necessary to go to higher order in . 
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3.3. Slow-roll, small-a regime 

We now reduce ourselves to the simplest case, in which we take both the slow-roll 
approximation and the small-a limit. Taking <C (f> <C 1, allows us to neglect all 
terms containing £2,3,4- The scalar field equation then implies that U' is also small, 
and so in the limit we are considering the dark energy must be mostly cosmological 
constant. Notice that to derive fl33l and flMj) we neglected terms of order a<p and, as a 
consequence, the correction turned out to be of order a 2 . It is therefore not correct to 
take the limit a — > in (1331) . (1341) and (1241) . (1251) have to be used instead. Equation (1241) 
in the slow-rolling regime gives 

7^ = l-«iil+ l^m) + 0(a 2 , 2 , 4>) . (36) 

Similarly from (l25j) we obtain, after some simplification, 

G, = G (1 + i±|^Q^ + 0(a 2 , 2 , 0) . (37) 

Differentiating and simplifying this gives a particularly simple relation 
C 91 

_^ = _ a m fi 1 + O(« 2 ,0 2 ,0). (38) 

Notice that both 7 and G*/G* are independent of the potential. The slow-roll and 
small-a limit allows us to put constraints directly on regardless of the potential U. 



4. Observing a coupled Gauss-Bonnet gravity 

We have seen that we can derive the two LPN observable quantities 7 and G*/G* in 
terms of the GB parameters. Now we discuss briefly how to use these quantities to 
put limits on the model. We will explicitly describe two methods: growth of matter 
perturbations and the integrated Sachs- Wolfe (ISW) effect. In both cases we will see that 
constraints from current data are rather weak, but non-trivial; moreover, the prospects 
for future experiments are encouraging. For simplicity, we consider only the limit of 
slow-rolling and small a in this section. 

The equation for the matter growth that we get from ( |371) is 

This is just the standard perturbation equation with an effective matter parameter 

fieff = U + On • (40) 

We can therefore obtain an approximate solution for the growth factor f = 5/5 in the 
standard way. In conventional gravity this is f — with the numerical coefficient 
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q ~ 0.6 [24J. To include the effects of the GB term, we simply replace fi m with fi eff . For 
small fii this implies 

i^ ai ( 1+q i±i^ a y (41) 



8 m V 2 

For typical values of fi m = 0.3 and q = 0.6, we see that the GB term increases the 
standard ACDM growth rate by the fraction « l.Sfii. Comparing with the observational 
result fo = 0.51 ±0.1 [25] we obtain for the present value —0.2 < fii < 0.4, as we show in 
Figured! There are many upcoming observational programmes to detect the evolution of 
clustering at various redshifts, so future data will no doubt strengthen these constraints 
and extend their temporal range. 

As a second observational constraint we consider the ISW effect, which depends on 
the variation of Newton's potential From ( 12111 we have for every wavenumber k 

I— 14 

where we used the relation Q m H 2 ~ p m a 2 ~ a l . Then we obtain 



p m a 
g(4 + 3fi m ) n 



y fi m fii • (43) 



2 

The standard ACDM result is then modified by a term that for fi m « 0.3 amounts 
to — 2.4fii. Comparing with the observational result is rather difficult here, because 
observational groups typically only quote their results in terms of the ISW cross- 
correlation with large scale structure. Hence knowledge would be required of the 
background dynamics, and of the full evolution of perturbations at all times, which 
is beyond the scope of this work. Moreover the uncertainties in the data so far are so 
large, even for ACDM models, that a very detailed comparison looks premature. A rough 
limit on Oi can be obtained taking the range of fi m values that are consistent with ISW 
measurements for a pure ACDM or for a constant-it? dark energy model. From [26 | l2~T l l2"8] 
we can conservatively summarise the allowed range as 0.05 < fi m < 0.5. We then see 
from Figure [U that fii has to be in the interval (—0.2,0.1) in order not to produce too 
large or too small a ISW cross- correlation. Apparently the constraints from the ISW 
effect are similar to those from the structure growth; however, let us remark that this 
procedure assumes implicitly that in our coupled GB model the evolution of large scale 
structure is similar to standard dark energy models, which is clearly unwarranted. 

Combining the two constraints and assuming fi m = 0.3 ±0.1 at the present, we see 
that current data put a weak but not totally irrelevant constraint on fii, of the order 
of | < 0.15. This shows, firstly that the coupled GB term is unlikely to provide the 
bulk of the dark energy, and secondly that future data have the potential to greatly 
improve this bound. Finally, the fact that fii has to be quite a bit smaller than unity 
is in agreement with our assumptions of small a and slow-rolling. Of course all of this 
says nothing about the value of fii in the past. 
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0.2 0.25 0.3 0.35 0.4 



Figure 1. In this plot we show the constraints on Sli from the ISW effect (within 
the dotted lines) and from the growth rate / of structures (grey region), assuming a 
present value Q TO — 0.3 ± 0.1. 



5. Conclusion 

We have investigated extensions of quintessence models in which part of the dark 
energy comes from quadratic-order gravity terms (and corresponding scalar kinetic 
terms). Even though the models contained no direct coupling of the quintessence scalar 
field to matter, the effective gravitational coupling (G>) in the model acquired a time 
dependence and the stress parameter 7 deviates from its general-relativistic value. We 
included all possible ghost-free second-order gravity terms in our model, since it is 
natural for them to all be present simultaneously. This can seen, for example, by 
considering a scalar field which arises from the compactification of an iV-torus. 

The time variation of G* and the value of 7 depend on the contribution of the 
higher-order gravity terms in the Friedmann equation, i.e. on the density fractions fij. 
This suggests that by taking them to be small, the constraints can be satisfied. If we 
take the coupling of the higher-order gravity terms (a) to be small, this is indeed the 
case. On the other hand if we suppose that the variation of the scalar field is small 
(which the fi« are all proportional to), then we find that 7 — 1 and G*/G* can still 
be large. This is due to the fact that the higher-order gravity terms can still have a 
significant effect on perturbations to the scalar field, even when they have little effect 
on the Friedmann equation. 

We considered the possibility of observing the coupled GB term through the matter 
perturbation growth and through the ISW effect. We found that current data limit 
I f2i| < 0.15, assuming both small a and slow rolling. This implies that the present 
effect of a coupled GB term is quite limited; on the other hand, its impact during other 
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epochs in cosmic history remains unbounded. 

The contributions from the higher-order "fc-essence" scalar kinetic terms (^3,4) 
only affect 7 and G* indirectly via the dynamics of the scalar field, and so the direct 
constraints on them are much weaker (and completely absent in our slow-roll limit). 
However, unless the couplings £3,4 are unnaturally large, the density fractions f2 3>4 will 
be smaller than fl^ by factors of <f) 2 ,<p 3 respectively, and thus completely negligible in 
a slowly rolling expansion. Moreover, we have seen that the slow-roll condition is the 
less restrictive one, in the sense that any faster dynamics will in general yield (barring 
chance cancellation) tighter constraints on (we have shown this explicitly only in the 
small-a limit, but we can conjecture it is a general property, since higher-order terms 
in a will generally also bring in higher-order terms in the time derivatives of <p). Of 
course it is possible that the corrections to G> coming from the different modifications 
to gravity will cancel each other, resulting in far weaker constraints. This will generally 
require fine-tuning, although it could occur naturally as a result of symmetries in the 
underlying theory. 

Finally, we emphasise that the constraints obtained here are restricted to the 
Einstein frame and to energy conservation. This hypothesis permitted us to disentangle 
any deviations from general relativity exclusively from the higher-order terms. Relaxing 
this hypothesis, and also the determination of exact cosmological solutions, are 
interesting subjects that we leave for future studj{JJ. 
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